Virus detection via programmable Type III-A CRISPR-Cas systems.

Among the currently available virus detection assays, those based on the programmable CRISPR-Cas enzymes have the advantage of rapid reporting and high sensitivity without the requirement of thermocyclers. Type III-A CRISPR-Cas system is a multi-component and multipronged immune effector, activated by viral RNA that previously has not been repurposed for disease detection owing in part to the complex enzyme reconstitution process and functionality. Here, we describe the construction and application of a virus detection method, based on an in vivo-reconstituted Type III-A CRISPR-Cas system. This system harnesses both RNA- and transcription-activated dual nucleic acid cleavage activities as well as internal signal amplification that allow virus detection with high sensitivity and at multiple settings. We demonstrate the use of the Type III-A system-based method in detection of SARS-CoV-2 that reached 2000 copies/µl sensitivity in amplification-free and 60 copies/µl sensitivity via isothermal amplification within 30 min and diagnosed SARS-CoV-2-infected patients in both settings. The high sensitivity, flexible reaction conditions, and the small molecular-driven amplification make the Type III-A system a potentially unique nucleic acid detection method with broad applications.

The native cistrome and sequence motif families of the maize ear.

Elucidating the transcriptional regulatory networks that underlie growth and development requires robust ways to define the complete set of transcription factor (TF) binding sites. Although TF-binding sites are known to be generally located within accessible chromatin regions (ACRs), pinpointing these DNA regulatory elements globally remains challenging. Current approaches primarily identify binding sites for a single TF (e.g. ChIP-seq), or globally detect ACRs but lack the resolution to consistently define TF-binding sites (e.g. DNAse-seq, ATAC-seq). To address this challenge, we developed MNase-defined cistrome-Occupancy Analysis (MOA-seq), a high-resolution (< 30 bp), high-throughput, and genome-wide strategy to globally identify putative TF-binding sites within ACRs. We used MOA-seq on developing maize ears as a proof of concept, able to define a cistrome of 145,000 MOA footprints (MFs). While a substantial majority (76%) of the known ATAC-seq ACRs intersected with the MFs, only a minority of MFs overlapped with the ATAC peaks, indicating that the majority of MFs were novel and not detected by ATAC-seq. MFs were associated with promoters and significantly enriched for TF-binding and long-range chromatin interaction sites, including for the well-characterized FASCIATED EAR4, KNOTTED1, and TEOSINTE BRANCHED1. Importantly, the MOA-seq strategy improved the spatial resolution of TF-binding prediction and allowed us to identify 215 motif families collectively distributed over more than 100,000 non-overlapping, putatively-occupied binding sites across the genome. Our study presents a simple, efficient, and high-resolution approach to identify putative TF footprints and binding motifs genome-wide, to ultimately define a native cistrome atlas.

Multiple roles of H2A.Z in regulating promoter chromatin architecture in human cells.

Chromatin accessibility of a promoter is fundamental in regulating transcriptional activity. The histone variant H2A.Z has been shown to contribute to this regulation, but its role has remained poorly understood. Here, we prepare high-depth maps of the position and accessibility of H2A.Z-containing nucleosomes for all human Pol II promoters in epithelial, mesenchymal and isogenic cancer cell lines. We find that, in contrast to the prevailing model, many different types of active and inactive promoter structures are observed that differ in their nucleosome organization and sensitivity to MNase digestion. Key aspects of an active chromatin structure include positioned H2A.Z MNase resistant nucleosomes upstream or downstream of the TSS, and a MNase sensitive nucleosome at the TSS. Furthermore, the loss of H2A.Z leads to a dramatic increase in the accessibility of transcription factor binding sites. Collectively, these results suggest that H2A.Z has multiple and distinct roles in regulating gene expression dependent upon its location in a promoter.

Chromatin structure profile data from DNS-seq: Differential nuclease sensitivity mapping of four reference tissues of B73 maize (Zea mays L).

Presented here are data from Next-Generation Sequencing of differential micrococcal nuclease digestions of formaldehyde-crosslinked chromatin in selected tissues of maize (Zea mays) inbred line B73. Supplemental materials include a wet-bench protocol for making DNS-seq libraries, the DNS-seq data processing pipeline for producing genome browser tracks. This report also includes the peak-calling pipeline using the iSeg algorithm to segment positive and negative peaks from the DNS-seq difference profiles. The data repository for the sequence data is the NCBI SRA, BioProject Accession PRJNA445708.

iSeg: an efficient algorithm for segmentation of genomic and epigenomic data.

BACKGROUND: Identification of functional elements of a genome often requires dividing a sequence of measurements along a genome into segments where adjacent segments have different properties, such as different mean values. Despite dozens of algorithms developed to address this problem in genomics research, methods with improved accuracy and speed are still needed to effectively tackle both existing and emerging genomic and epigenomic segmentation problems. RESULTS: We designed an efficient algorithm, called iSeg, for segmentation of genomic and epigenomic profiles. iSeg first utilizes dynamic programming to identify candidate segments and test for significance. It then uses a novel data structure based on two coupled balanced binary trees to detect overlapping significant segments and update them simultaneously during searching and refinement stages. Refinement and merging of significant segments are performed at the end to generate the final set of segments. By using an objective function based on the p-values of the segments, the algorithm can serve as a general computational framework to be combined with different assumptions on the distributions of the data. As a general segmentation method, it can segment different types of genomic and epigenomic data, such as DNA copy number variation, nucleosome occupancy, nuclease sensitivity, and differential nuclease sensitivity data. Using simple t-tests to compute p-values across multiple datasets of different types, we evaluate iSeg using both simulated and experimental datasets and show that it performs satisfactorily when compared with some other popular methods, which often employ more sophisticated statistical models. Implemented in C++, iSeg is also very computationally efficient, well suited for large numbers of input profiles and data with very long sequences. CONCLUSIONS: We have developed an efficient general-purpose segmentation tool and showed that it had comparable or more accurate results than many of the most popular segment-calling algorithms used in contemporary genomic data analysis. iSeg is capable of analyzing datasets that have both positive and negative values. Tunable parameters allow users to readily adjust the statistical stringency to best match the biological nature of individual datasets, including widely or sparsely mapped genomic datasets or those with non-normal distributions.

Regulated large-scale nucleosome density patterns and precise nucleosome positioning correlate with V(D)J recombination.

We show that the physical distribution of nucleosomes at antigen receptor loci is subject to regulated cell type-specific and lineage-specific positioning and correlates with the accessibility of these gene segments to recombination. At the Ig heavy chain locus (IgH), a nucleosome in pro-B cells is generally positioned over each IgH variable (VH) coding segment, directly adjacent to the recombination signal sequence (RSS), placing the RSS in a position accessible to the recombination activating gene (RAG) recombinase. These changes result in establishment of a specific chromatin organization at the RSS that facilitates accessibility of the genomic DNA for the RAG recombinase. In contrast, in mouse embryonic fibroblasts the coding segment is depleted of nucleosomes, which instead cover the RSS, thereby rendering it inaccessible. Pro-T cells exhibit a pattern intermediate between pro-B cells and mouse embryonic fibroblasts. We also find large-scale variations of nucleosome density over hundreds of kilobases, delineating chromosomal domains within IgH, in a cell type-dependent manner. These findings suggest that developmentally regulated changes in nucleosome location and occupancy, in addition to the known chromatin modifications, play a fundamental role in regulating V(D)J recombination. Nucleosome positioning-which has previously been observed to vary locally at individual enhancers and promoters-may be a more general mechanism by which cells can regulate the accessibility of the genome during development, at scales ranging from several hundred base pairs to many kilobases.

Label-Free Relative Quantitation of Isobaric and Isomeric Human Histone H2A and H2B Variants by Fourier Transform Ion Cyclotron Resonance Top-Down MS/MS.

Histone variants are known to play a central role in genome regulation and maintenance. However, many variants are inaccessible by antibody-based methods or bottom-up tandem mass spectrometry due to their highly similar sequences. For many, the only tractable approach is with intact protein top-down tandem mass spectrometry. Here, ultra-high-resolution FT-ICR MS and MS/MS yield quantitative relative abundances of all detected HeLa H2A and H2B isobaric and isomeric variants with a label-free approach. We extend the analysis to identify and relatively quantitate 16 proteoforms from 12 sequence variants of histone H2A and 10 proteoforms of histone H2B from three other cell lines: human embryonic stem cells (WA09), U937, and a prostate cancer cell line LaZ. The top-down MS/MS approach provides a path forward for more extensive elucidation of the biological role of many previously unstudied histone variants and post-translational modifications.

Hierarchical regulation of the genome: global changes in nucleosome organization potentiate genome response.

Nucleosome occupancy is critically important in regulating access to the eukaryotic genome. Few studies in human cells have measured genome-wide nucleosome distributions at high temporal resolution during a response to a common stimulus. We measured nucleosome distributions at high temporal resolution following Kaposi's-sarcoma-associated herpesvirus (KSHV) reactivation using our newly developed mTSS-seq technology, which maps nucleosome distribution at the transcription start sites (TSS) of all human genes. Nucleosomes underwent widespread changes in organization 24 hours after KSHV reactivation and returned to their basal nucleosomal architecture 48 hours after KSHV reactivation. The widespread changes consisted of an indiscriminate remodeling event resulting in the loss of nucleosome rotational phasing signals. Additionally, one in six TSSs in the human genome possessed nucleosomes that are translationally remodeled. 72% of the loci with translationally remodeled nucleosomes have nucleosomes that moved to positions encoded by the underlying DNA sequence. Finally we demonstrated that these widespread alterations in nucleosomal architecture potentiated regulatory factor binding. These descriptions of nucleosomal architecture changes provide a new framework for understanding the role of chromatin in the genomic response, and have allowed us to propose a hierarchical model for chromatin-based regulation of genome response.

Comprehensive nucleosome mapping of the human genome in cancer progression.

Altered chromatin structure is a hallmark of cancer, and inappropriate regulation of chromatin structure may represent the origin of transformation. Important studies have mapped human nucleosome distributions genome wide, but the role of chromatin structure in cancer progression has not been addressed. We developed a MNase-Transcription Start Site Sequence Capture method (mTSS-seq) to map the nucleosome distribution at human transcription start sites genome-wide in primary human lung and colon adenocarcinoma tissue. Here, we confirm that nucleosome redistribution is an early, widespread event in lung (LAC) and colon (CRC) adenocarcinoma. These altered nucleosome architectures are consistent between LAC and CRC patient samples indicating that they may serve as important early adenocarcinoma markers. We demonstrate that the nucleosome alterations are driven by the underlying DNA sequence and potentiate transcription factor binding. We conclude that DNA-directed nucleosome redistributions are widespread early in cancer progression. We have proposed an entirely new hierarchical model for chromatin-mediated genome regulation.

Stimulation of the Drosophila immune system alters genome-wide nucleosome occupancy.

In eukaryotes, nucleosomes participate in all DNA-templated events by regulating access to the underlying DNA sequence. However, nucleosome dynamics during a genome response have not been well characterized [1,2]. We stimulated Drosophila S2 cells with heat-killed Gram-negative bacteria Salmonella typhimurium, and mapped genome-wide nucleosome occupancy at high temporal resolution by MNase-seq using Illumina HiSeq 2500. We show widespread nucleosome occupancy change in S2 cells during the immune response, with the significant nucleosomal loss occurring at 4 h after stimulation. Data have been deposited to the Gene Expression Omnibus (GEO) database repository with the dataset identifier GSE64507.

Nucleosome Repositioning: A Novel Mechanism for Nicotine- and Cocaine-Induced Epigenetic Changes.

Drugs of abuse modify behavior by altering gene expression in the brain. Gene expression can be regulated by changes in DNA methylation as well as by histone modifications, which alter chromatin structure, DNA compaction and DNA accessibility. In order to better understand the molecular mechanisms directing drug-induced changes in chromatin structure, we examined DNA-nucleosome interactions within promoter regions of 858 genes in human neuroblastoma cells (SH-SY5Y) exposed to nicotine or cocaine. Widespread, drug- and time-resolved repositioning of nucleosomes was identified at the transcription start site and promoter region of multiple genes. Nicotine and cocaine produced unique and shared changes in terms of the numbers and types of genes affected, as well as repositioning of nucleosomes at sites which could increase or decrease the probability of gene expression based on DNA accessibility. Half of the drug-induced nucleosome positions approximated a theoretical model of nucleosome occupancy based on physical and chemical characteristics of the DNA sequence, whereas the basal or drug naïve positions were generally DNA sequence independent. Thus we suggest that nucleosome repositioning represents an initial dynamic genome-wide alteration of the transcriptional landscape preceding more selective downstream transcriptional reprogramming, which ultimately characterizes the cell- and tissue-specific responses to drugs of abuse.

BAF250a Protein Regulates Nucleosome Occupancy and Histone Modifications in Priming Embryonic Stem Cell Differentiation.

The unique chromatin signature of ES cells is fundamental to the pluripotency and differentiation of ES cells. One key feature is the poised chromatin state of master developmental genes that are transcriptionally repressed in ES cells but ready to be activated in response to differentiation signals. Poised chromatin in ES cells contains both H3 Lys-4 trimethylation (H3K4me3) and H3 Lys-27 trimethylation (H3K27me3) methylation, indicating activating and repressing potential. However, the contribution of non-covalent chromatin structure to the poised state is not well understood. To address whether remodeling of nucleosomes is important to the poised state, we characterized the function of BAF250a, a key regulatory subunit of the ES cell ATP-dependent Brahma-associated factor (BAF) chromatin remodeling complex (esBAF). Acute deletion of BAF250a disrupted the differentiation potential of ES cells by altering the expression timing of key developmental genes and pluripotent genes. Our genome-wide nucleosome and histone modification analyses indicated that the disruption of gene expression timing was largely due to changes of chromatin structures at poised genes, particularly those key developmental genes mediated by BAF250a. Specifically, BAF250a deletion caused a nucleosome occupancy increase at H3K4me3- and/or H3K27me3-associated promoters. Moreover, H3K27me3 levels and the number of bivalent promoter genes were reduced in BAF250a KO ES cells. We revealed that BAF250a ablation led to elevated Brg1 but reduced Suz12 recruitment at nucleosome occupancy-increased regions, indicating an unexpected and complicated role of BAF250a in regulating esBAF and Polycomb repressive complex (PRC) activities. Together, our studies identified that BAF250a mediates esBAF and PRC functions to establish the poised chromatin configuration in ES cells, which is essential for the proper differentiation of ES cells.

DNA-Encoded Chromatin Structural Intron Boundary Signals Identify Conserved Genes with Common Function.

The regulation of metazoan gene expression occurs in part by pre-mRNA splicing into mature RNAs. Signals affecting the efficiency and specificity with which introns are removed have not been completely elucidated. Splicing likely occurs cotranscriptionally, with chromatin structure playing a key regulatory role. We calculated DNA encoded nucleosome occupancy likelihood (NOL) scores at the boundaries between introns and exons across five metazoan species. We found that (i) NOL scores reveal a sequence-based feature at the introns on both sides of the intron-exon boundary; (ii) this feature is not part of any recognizable consensus sequence; (iii) this feature is conserved throughout metazoa; (iv) this feature is enriched in genes sharing similar functions: ATPase activity, ATP binding, helicase activity, and motor activity; (v) genes with these functions exhibit different genomic characteristics; (vi) in vivo nucleosome positioning data confirm ontological enrichment at this feature; and (vii) genes with this feature exhibit unique dinucleotide distributions at the intron-exon boundary. The NOL scores point toward a physical property of DNA that may play a role in the mechanism of pre-mRNA splicing. These results provide a foundation for identification of a new set of regulatory DNA elements involved in splicing regulation.

Topologically associating domains are stable units of replication-timing regulation.

Eukaryotic chromosomes replicate in a temporal order known as the replication-timing program. In mammals, replication timing is cell-type-specific with at least half the genome switching replication timing during development, primarily in units of 400-800 kilobases ('replication domains'), whose positions are preserved in different cell types, conserved between species, and appear to confine long-range effects of chromosome rearrangements. Early and late replication correlate, respectively, with open and closed three-dimensional chromatin compartments identified by high-resolution chromosome conformation capture (Hi-C), and, to a lesser extent, late replication correlates with lamina-associated domains (LADs). Recent Hi-C mapping has unveiled substructure within chromatin compartments called topologically associating domains (TADs) that are largely conserved in their positions between cell types and are similar in size to replication domains. However, TADs can be further sub-stratified into smaller domains, challenging the significance of structures at any particular scale. Moreover, attempts to reconcile TADs and LADs to replication-timing data have not revealed a common, underlying domain structure. Here we localize boundaries of replication domains to the early-replicating border of replication-timing transitions and map their positions in 18 human and 13 mouse cell types. We demonstrate that, collectively, replication domain boundaries share a near one-to-one correlation with TAD boundaries, whereas within a cell type, adjacent TADs that replicate at similar times obscure replication domain boundaries, largely accounting for the previously reported lack of alignment. Moreover, cell-type-specific replication timing of TADs partitions the genome into two large-scale sub-nuclear compartments revealing that replication-timing transitions are indistinguishable from late-replicating regions in chromatin composition and lamina association and accounting for the reduced correlation of replication timing to LADs and heterochromatin. Our results reconcile cell-type-specific sub-nuclear compartmentalization and replication timing with developmentally stable structural domains and offer a unified model for large-scale chromosome structure and function.

Differential nuclease sensitivity profiling of chromatin reveals biochemical footprints coupled to gene expression and functional DNA elements in maize.

The eukaryotic genome is organized into nucleosomes, the fundamental units of chromatin. The positions of nucleosomes on DNA regulate protein-DNA interactions and in turn influence DNA-templated events. Despite the increasing number of genome-wide maps of nucleosome position, how global changes in gene expression relate to changes in nucleosome position is poorly understood. We show that in nucleosome occupancy mapping experiments in maize (Zea mays), particular genomic regions are highly susceptible to variation introduced by differences in the extent to which chromatin is digested with micrococcal nuclease (MNase). We exploited this digestion-linked variation to identify protein footprints that are hypersensitive to MNase digestion, an approach we term differential nuclease sensitivity profiling (DNS-chip). Hypersensitive footprints were enriched at the 5' and 3' ends of genes, associated with gene expression levels, and significantly overlapped with conserved noncoding sequences and the binding sites of the transcription factor KNOTTED1. We also found that the tissue-specific regulation of gene expression was linked to tissue-specific hypersensitive footprints. These results reveal biochemical features of nucleosome organization that correlate with gene expression levels and colocalize with functional DNA elements. This approach to chromatin profiling should be broadly applicable to other species and should shed light on the relationships among chromatin organization, protein-DNA interactions, and genome regulation.

Changes in nucleosome occupancy occur in a chromosome specific manner.

In the eukaryotic nucleus, DNA is packaged into chromatin. The fundamental subunit of chromatin is the nucleosome, DNA is wrapped 1.6 times around a histone octamer core. Nuclear processes in eukaryotes are impacted by whether regulatory DNA is occupied by nucleosomes. We used microarrays to measure nucleosome occupancy in human cells post Kaposi's sarcoma-associated herpesvirus (KSHV) reactivation at hundreds of immunity-related loci. The detailed analysis of these technologies can be found in recent publications from our lab (Druliner et al., 2013; Sexton et al., 2014). We found that nucleosome redistributions displayed chromosome specific nucleosome occupancy. This resource can be used to map nucleosome distributions in a variety of biological contexts.

The spring-loaded genome: nucleosome redistributions are widespread, transient, and DNA-directed.

Nucleosome occupancy plays a key role in regulating access to eukaryotic genomes. Although various chromatin regulatory complexes are known to regulate nucleosome occupancy, the role of DNA sequence in this regulation remains unclear, particularly in mammals. To address this problem, we measured nucleosome distribution at high temporal resolution in human cells at hundreds of genes during the reactivation of Kaposi's sarcoma-associated herpesvirus (KSHV). We show that nucleosome redistribution peaks at 24 h post-KSHV reactivation and that the nucleosomal redistributions are widespread and transient. To clarify the role of DNA sequence in these nucleosomal redistributions, we compared the genes with altered nucleosome distribution to a sequence-based computer model and in vitro-assembled nucleosomes. We demonstrate that both the predicted model and the assembled nucleosome distributions are concordant with the majority of nucleosome redistributions at 24 h post-KSHV reactivation. We suggest a model in which loci are held in an unfavorable chromatin architecture and "spring" to a transient intermediate state directed by DNA sequence information. We propose that DNA sequence plays a more considerable role in the regulation of nucleosome positions than was previously appreciated. The surprising findings that nucleosome redistributions are widespread, transient, and DNA-directed shift the current perspective regarding regulation of nucleosome distribution in humans.

Chromatin patterns associated with lung adenocarcinoma progression.

The development and progression of lung adenocarcinoma, one of the most common cancers, is driven by the interplay of genetic and epigenetic changes and the role of chromatin structure in malignant transformation remains poorly understood. We used systematic nucleosome distribution and chromatin accessibility microarray mapping platforms to analyze the genome-wide chromatin structure from normal tissues and from primary lung adenocarcinoma of different grades and stages. We identified chromatin-based patterns across different patients with lung adenocarcinoma of different cancer grade and stage. Low-grade cancers had nucleosome distributions very different compared with the corresponding normal tissue but had nearly identical chromatin accessibility. Conversely, nucleosome distributions of high-grade cancers showed few differences. Substantial disruptions in chromosomal accessibility were seen in a patient with a high-grade and high-stage tumor. These data imply that chromatin structure changes during the progression of lung adenocarcinoma. We have therefore developed a model in which low-grade lung adenocarcinomas are linked to changes in nucleosome distributions, whereas higher-grade tumors are linked to large-scale chromosomal changes. These results provide a foundation for the development of a comprehensive framework linking the general and locus-specific roles of chromatin structure to lung cancer progression. We propose that this strategy has the potential to identify a new class of chromatin-based diagnostic, prognostic and therapeutic markers in cancer progression.

Genome-wide prediction of nucleosome occupancy in maize reveals plant chromatin structural features at genes and other elements at multiple scales.

The nucleosome is a fundamental structural and functional chromatin unit that affects nearly all DNA-templated events in eukaryotic genomes. It is also a biochemical substrate for higher order, cis-acting gene expression codes and the monomeric structural unit for chromatin packaging at multiple scales. To predict the nucleosome landscape of a model plant genome, we used a support vector machine computational algorithm trained on human chromatin to predict the nucleosome occupancy likelihood (NOL) across the maize (Zea mays) genome. Experimentally validated NOL plots provide a novel genomic annotation that highlights gene structures, repetitive elements, and chromosome-scale domains likely to reflect regional gene density. We established a new genome browser (http://www.genomaize.org) for viewing support vector machine-based NOL scores. This annotation provides sequence-based comprehensive coverage across the entire genome, including repetitive genomic regions typically excluded from experimental genomics data. We find that transposable elements often displayed family-specific NOL profiles that included distinct regions, especially near their termini, predicted to have strong affinities for nucleosomes. We examined transcription start site consensus NOL plots for maize gene sets and discovered that most maize genes display a typical +1 nucleosome positioning signal just downstream of the start site but not upstream. This overall lack of a -1 nucleosome positioning signal was also predicted by our method for Arabidopsis (Arabidopsis thaliana) genes and verified by additional analysis of previously published Arabidopsis MNase-Seq data, revealing a general feature of plant promoters. Our study advances plant chromatin research by defining the potential contribution of the DNA sequence to observed nucleosome positioning and provides an invariant baseline annotation against which other genomic data can be compared.